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ABSTRACT 


A computer data reduction system has been reprogrammed to 
sample individual scan recordings of the line-spread function 
of laser light propagating through the turbulent marine boud- 
ary layer. 

The distribution of spatial frequencies of the Modulation 
Transfer Functions from random single scan samples from 
experimental data sets has been investigated. A preliminary 
analysis indicates that the zero spatial frequency component 


of the MIF is log-normal distributed. 
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I. INTRODUCTION 


Atmospheric turbulence adversely affects the propagation 
of a laser beam or the performance of an imaging system. The 
Modulation Transfer Function (MIF) is a measure of the degra- 
dation of resolution and thus the relative ability of the 
atmosphere to transmit spatial frequency information. Knowl- 
edge of the MIF of the atmosphere is necessary for prediction 
of the beam cross-sectional profile on target and spot wander 
for a high energy laser and its beam projection system. 

Members of the faculty of the Naval Postgraduate School 
(NPS), Department of Physics and Chemistry, have investigated 
the effects of atmospheric turbulence on the propagation of 
optical radiation in the marine boundary layer. 

The work started in 1973 and has continued through 1978. 
The research has been concentrated on developing a model for 
the MTF depending on wavelength, range and a properly path- 
position weighted value of the turbulence structure constant 
for optical index Cn. The work for the period 1 July 1974 - 
30 September 1977 is described in Ref. 1. Field experiments 
were performed by transmitting laser light of different wave- 
lengths along various seapaths. The received intensity was 
detected by a slit-scanning telescope system and recorded on 
tape. The Cn value was determined by computer calculations. 

The objective of this report is to describe the develop- 


ment of a computer program which samples individual scan 








recordings of the line-spread function and determines a pre- 
liminary statistical distribution of the spatial frequency 


content of the Modulation Transfer Function. 


LO 








IT. ANALYSIS OF DATA 


A. FUNCTIONAL DESCRIPTION 
Ll. Data 

Experimental data are the intensities of the line- 
spread functions as recorded on a 4-channel magnetic tape. 
One track is used for verbal identification of date, time 
and pertinent experimental information. Two tracks are used 
for recording scans and the last track is used for recording 
the trigger signal. 

Verbal information has been cross-checked against 
experiment log book and MTF data output. Data have been 
recorded for one and two wavelengths simultaneously. Only 
data for .6328 microns have been analyzed in this report. 

The line-spread function was sampled by the scanning detector 
at a rate of 50Hz and recorded at a tape speed of 37.5 inches 
per second. 

2. Data Reduction 

The data reduction is done on a Nicolet Instrument 
Company, NIC-80 data processor. This is a data digitizer 
followed by a 12K 20-bit digital computer and a Tektronix CRT 
display unit. The NIC-80 computer is controlled by a Hewlett- 
Packard HP-9825 desk top type programmable calculator. 

The NIC-80 computer reads the line-spread function 
for one scan, performs a Fourier Transform (FT) analysis of 
the signal and transfers the FT spatial frequency components 


O - 19 to the HP-9825. The MTF of the atmosphere is calculated 


ies 





by dividing the FT of the line-spread function by the FT of 
the optics. The FT of the optics is stored by the NIC-80 and 
transfered to HP-9825. The experimental analog signal is 
digitized at 1024 points. The DC component is removed by 
manual setting of signal bias based upon observation of scan 
on CRT display. Any remaining DC component is compensated by 
subtracting the average intensity value of 32 points at each 
end of the digitized signal. 

When not performing FT and transfer operations, the 
NIC-80 remains in signal display mode with successive scans 
being displayed on CRT. 

The computation and transfer time for one scan has 
been observed to be 5 seconds, and will be the time between 
sampling of scans. During data reduction, the tape recorder 
was operated at 3.75 IPS, corresponding to 5 scans per second. 
This implies that about every 25th scan of the recorded data 
was sampled for analysis. The number of scans sampled was 
limited by the length of experimental data recorded, normally 
in the order of one minute, or 3000 scans. This allowed a 
sample size of 100 - 125 for data analysis. An attempt to 
run the recorder at .375 IPS speed was not successful due to 
problems with the trigger signal. 

3. Statistical Analysis 

The HP-9825 calculates the MIF at 20 different spatial 
frequencies from 0 to 38000 per radian in steps of 2000 per 
radian. These spatial frequency points are labeled 0 to 19 


and are tested for normal and log-normal distribution. Goodness 


IZ 








of fit is determined by the Chi Square test. A data summary 
containing mean, standard deviation and Chi Square values 
for normal and log-normal tests is outputted on the printer. 
Mean and + 1 standard deviation for frequency points 0 - 19 
is plotted, with the zero spatial frequency normalized to a 
value of unity and all other data adjusted correspondingly. 
A histogram showing relative distribution of MTF values for 
one specific spatial frequency is plotted as requested in 
program input. 

Examples of the above are given in Figures 1, 2 and 3 
respectively. The relationship of the MTF distribution to 
the Cn value has been investigated. The Cn values ranged in 


: to Lxl07/. Typical line-spread functions 


the order of 8x10 
for two different Cn values are shown in Figure 4. Note that 
Cn squared values are given. 

The value of Cn used has been taken from existing MTF 
data outputs. If no output was available, an MIF analysis 
was performed. A detailed discussion of MIF and Cn is given 
in Ref. 1. The scan presentation is a function of time and 
detector position. The dwell time setting was obtained by 
analyzing calibration data recorded on tape. Dwell time is 
the sampling time per data point. Data have been obtained 
by recording the detector output with a diffraction grating 
of known dimensions placed in front of the telescope. This 
provides a known reference of angular spread. The dwell time 
ls adjusted so that the angular spread of the first order 
diffraction as observed on the display is equal to the ref- 


erence value. 
13 
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B. PROGRAM THEORETICAL CONSIDERATIONS 
The data were tested for normal and in normal distribucion. 
L. Normal Curve 


The normal or Gaussian curve is expressed as 


Y = constant exp (-& (A)? 


) 


where 


x is the value of MIF for normal test distribution 
or value of In MTF for log-normal distribution. 


x 1S fhe Uunpiased estimate of the Mean of the MIF 
values or Mean of ln MTF values. 


S is the unbiased estimate of the standard deviation 
Ot Phe dL serEipDUuLion. 


Y is the probability density of the distribution. 


Constant = i/svY27, where i is the width of a histogram 
bin to make area under curve equal to li. 


The Mean x is the arithmetic Mean given by: 
X = 


where N is the number of observations (samples). 
The standard deviation s is the positive root of the 


Variance s squared given by 


a 
S — 


lu 


(x; -&)°/ (N-1) 
a 


For computational purposes the following algebraically 
equivalent formula has been used: 
) N N 


s© = (5 (x,)*-( 5 x,)°/N)/(N-1) 
i=l °* i=l * 
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The probability density curve has been computed and 

plotted over a range of Mean + 3 standard deviations. 
2. Goodness of Fit 

Test of goodness of fit to tested distribution has 
been done by histogram and Chi Square calculation. 

The width of the histogram bins is a program variable 
expressed as a ratio of the standard deviation. 1/3 has been 
used corresponding to a total of 18 bins. The number of 
samples with values within a bin is expressed as a fraction 
of 1, so that the area of the histogram is equal to the width 
of the bin interval. 

The Chi Square goodness of fit value has been calcu- 


vated from 


2 
_ 5 (NE(n)-F(n)] 
=< ” “—NE(hny 
i=l 
where 
N = number of observations 
n = number of histogram bins 
F(n) = number of observations within a bin 


The term Nf(n) is the theoretical number of observations 
which should fall within a bin. The probability area has been 
calculated as the sum of a rectangle and a triangle approxi- 
mating the area under the normal curve. The sampling 
distribution of the Chi Square statistics is tabulated as a 
function of degrees of freedom and percentile. The number of 
degrees of freedom is calculated from 


dé = k-n-1 
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where 


k 


number of bins in histogram 

n = number of parameters estimated 
For a 18 bin histogram with mean and standard deviation 
estimated, the degrees of freedom is 15, with a 95 percentile 
value tabulated as 25. A calculated Chi Square value less 
than 25 will, at the 5% significance level, imply that the 


data fits the given distribution. 


C. COMPUTER PROGRAMS 
l. NIC-80 Data Reduction Program 
This program has been developed utilizing subroutines 
from the existing MIF analysis program. Pages 3 and 4 of the 
MIF program have not been changed. Main program is on page l 
and a data transfer subroutine has been inserted on page 2. 
Program flowchart is given in Figure 5A. Program listing is 
ene pace LZ. 
2. Data Analysis Program 
Due to limited HP-9825 memory capacity, this program 
has been split into two parts stored on separate program 
files. During operation, programs are loaded and executed 
under program control. Program ‘'A'' is on file 8 and contains 
NIC-80 control and all analysis functions except data summary, 
which is on program "B" located on file 9. Programs are 
documented on flowcharts (Figures 6 and 7). Program listings 
are on pages 117-121. Execution of analysis program requires 


HP tape inserted with track 1 available for data storage. 


20 





Each data file must be minimum 160 bytes long and the number 
of files must be equal to or larger than the number of scans 
sampled. 

Execution time for a 100 scan analysis run is of the 
order of 90 minutes. Most of this time is required for tape 


movements and operations. 
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Figure 5A. Program Organization 
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Figure 5B. Program Organization 
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Figure 6A. Data Analysis Program 
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Figure 6B. Data Analysis Program 
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Figure 6D. Data Analysis Program 
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Figure 7A. Data Summary Program 
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Figure 8B. Data Summary Program 
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Lift. RESULIS 


A total of 24 data analysis runs have been performed. 
Five of these have been on the same data sets in order to 
investigate consistency of distribution within the same data 
Bec . 

An arbitrary grading system for goodness of fit over a 
range of spatial frequencies has been established as follows: 
- All components in range within 95 percentile = l. 
~ One component in range outside 95 percentile = 2. 


~ Two or more components in range outside 95 percentile 


Data summary and MIF mean and standard deviation plots 
are given on pages 39-86, arranged in order of decreasing 
order of Cn. A distribution analysis summary is tabulated 
on pages 31-33. A set of histograms for all spatial frequen- 


cies for one data set is given on pages 87-112. 


A. CONSISTENCY OF DISTRIBUTION 
Data set number 12 has been run 5 times to verify consis- 
tency of distribution as only every 25th scan was sampled. 
Utilizing the grading system for spatial frequency ranges 
0, 1-6, 7-12 and 13-19, the following results as given in 


Table II and III, were found: 
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TABLE IL 
DISTIRUBTION ANALYSIS SUMMARY 


: Spat. Log 













1.04] 031176 | 0 0 freq. log 
| LG data best fit 
yy | ey) 
10 | Wee 19 







O freq. normal 
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deviation 


1.061 021877. | 1225 0 
ier 

x | jee 
= 13500 
1.06) 040276 | 1420 0 
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7 ao) 
LO } 13-19 
1.07] 031176 | 2247 0 
16 

a7 7-12 
10 13219 
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_7 Fano) 
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PG 
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Legend: Ll. All values within Chi Square 95 percentile. 
2. One value above Chi Square 95 percentile. 
3. Two or more values above Chi Square 
95 percentile. 


Good fit at 
low frequen- 
cles 







O,2neq... 16s 
data best fit 





O freq. log 
data best fit 


Baas Hetens tare bowen lores | eee 





rr’Wr lo WW - LO Po PF wi Lo Lo LO fo LO Wo rr ot? r re 


log data best 
i 


UO NO LO GO WI fo 


eee 
WI Lo 








DISTRIBUTION ANALYSIS SUMMARY 


cn Normown | Comments 
5.03| 092377 | 1452 0 
_ 
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10 13-19 
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-9 7-12 
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1-6 
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6.08 | 092277 | 1705 0 
1-6 
-9 7-12 
10 13-19 
6.92 | 121776 | 1012 0 
ee 
-8 7-12 
10 13-19 
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Gia99 021577 1700 
1-6 
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value 26.6 


! 
co 


log data low 
chi sq. values 


WNe Ee Lo oO NF 


Absolute values 
low 


1 This data set 

al run 5 times to 

1 check consis- 

Ll tency within 
data set 


Lo LO FR fo ae NO Go F bo 


Wn 
co 
LO 


Lo Lo Lo Lo GO Lo bo LO PO DY NOS 


aS eH NON Rr Lo ho OS fo 





log data range 
0-8 low chi sq. 
alues 








log data good 
5a 


WWE [Ree rE 


Cane 
a aan 


-g 7-12 
10 13-19 


Legend: 1. All values within Chi Square 95 percentile. 
2. One value above Chi Square 95 percentile. 
3. Two or more values above Chi Square 
95 percentile. 
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Legend: Il. All values within Chi Square 95 percentile. 
2. One value above Chi Square 95 percentile. 
3. Two or more values above Chi Square 
95 percentile. 
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TABLE II 


Grading Consistency of Data from Same Data Set 


Number of Runs with Grading 


Frequency Range Normal Log Normal 
Grade |e aes i ae: 
0 L 4 Q 4 1 0 
Ll - 6 1 0 4 See. 6 ame 8, 
7 - 12 CO a> i ees Rae 
13 - 19 ONO Ue cee ee 8. 


The number of samples which met the 95 percentile criterion 
for each spatial frequency in the O - 6 frequency range is 


fisted in Table IIL. 


TABLE II1L 


Consistency of Individual Frequency Components 
from Samples of Same Data Set 


Number of Samples within 95% 


Frequency Normal Log Normal 

0 il 4 

1 5 
2 2 5 
3 3 » 
4 3 5 
5 Z 5 
6 3 5 
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The results indicate that the distribution can be expected to 
give same Chi Square values for different sample series from 


same data set. 


B. FACTORS AFFECTING DISTRIBUTIONS 
im, ocimtrllaciron 
If the intensity of a laser beam which has propagated 
through a turbulent medium is observed at a point, the inten- 
sity appears to fluctuate or scintillate in a random way. 
It has been shown that the values of the logarithm 
of the amplitude follow a Gaussian distribution (Ref. l, 
page 15). Spatial frequency component zero represents the 
area under the curve (signal). This area is proportional to 
the received intensity. It is thus expected that the distri- 
bution of the zero spatial frequency component should fit 
tne log normal distribution. 
2. Signal to Noise Ratio 
The amplitude of the spatial frequencies decreases 
with increasing frequency number. For spatial frequency 15 
and above, the amplitude is of the order of 10% of the zero 
spatial frequency component. At higher frequencies noise 
will contribute more relative to the amplitude. Noise is 
considered normal distributed, and it is expected that the 


distribution at higher frequencies is normal. 


C. RESULTS AND DISCUSSION 
1. Log-Normal Distribution of 0 Spatial Frequency 
In 13 of 20 runs the spatial frequency zero fits the 
In normal Chi Square 95 percentile criterion. For runs which 
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do not meet the distribution criterion, the ratio of mean to 
standard deviation is in the order of 2. This may indicate 
saturation, and thus not follow log normal distribution (Ref. 
1, page 16). 
The results do indicate that the zero spatial freq- 
Hency is log normal distributed. 
2. Normal Distribution of Higher Frequencies 
For data analysis the modulus, e.g., absolute value 
of the data has been used, and sign of the data has not been 
included. This will introduce a skewness of distribution. 
It is believed that this is the reason for lack of normal 
distribution for data at higher frequencies. 
3. Distribution versus Cn Value 
Data have been arbitrarily grouped into ranges with 
respect to Cn values as follows: 


Cn High: Cn > 6210s" 


8 8 


Cn Medium: Cn between 2.10 ° and 6.10 


Ch Low: Cn < 2.1078 
The vatio of distributions graded 1 for the various 
Cn ranges is listed in Table IV. 
The results do not indicate that the distribution 
does change with Cn values. For all analysis results, it 


should be noted that the sample size is small, and the results 


should be interpreted as indications only. 
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TABLE IV 


Ratio of Distribution Graded 
1 for Different Cn Ranges 


Spacral 
Cn Frequency Normal Log Normal 
H 0 5/6 5/6 
I Ley 2/6 3/6 
G 7 = 12 3/6 2/6 
H 13 - 19 0/6 0/6 
M 0 S/T 4/7 
. lL - 6 ys S/d 
; 7 - 12 0/7 3/7 
‘ 13 - 19 0/7 Zi 7 
0 4/7 4/7 
1 - 6 Si 7 a7 
; (uo aZ 0/7 1/7 
. 13-2 17 Zi 7 


D. RECOMMENDATIONS FOR FURTHER WORK 
To reduce analysis computing time, the HP-9825 should be 
equipped with additional memory capacity. To eliminate use 
of tape for data storage, memory should be increased to 25K. 
The data analysis for each run should be stored and 
analyzed for mean and standard deviation. Several runs on 


the same data set should be performed to verify distribution. 


a7 





Number of scans sampled in each run should be increased 
by either solving signal trigger problem at low tape speed 
or by reprogramming NIC-80 computer to perform and store 


partial analysis results. 
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APPENDIX A 


Operating Instructions for Data Analysis 


A. Initialization 
1. NIC-80, Analysis program not loaded 

(a) Turn on computer with key, depress "STOP" 

(b) Turn on Data Bus unit 

(c) Set toggle switches 0 - Il to up position 

(d) Depress "LOAD PC" 

(e) Depress "CONTINUE" 

(f) Press "EXECUTE" 

2. HP-9825 

(a) Turn on HP-9825 

(b) Insert Data Analysis cassette and load track 0, 
file 1 for NIC-80 loading program 

Ce) Push RUN? = 3c, di, dg bet, Oe Sn God Neo laved. 

(d) Push "n" and "CONTINUE"; load file # is displayed. 
The NIC-80 data analysis program is located on files 2, 3, 4 
and 5. Enter one of these numbers and push "CONTINUE". 

Ce) Display Ce Ga is ae ee Oe Ge ee PUG 
and continue. The section of the program will be transferred 
to NIC-80. At the end of each transfer, check sums are 
printed and should be identical. 

(f) Repeat steps (d) and (e) until entire program is 
transferred. If NIC-80 contained MIF analysis program pre- 


viously, only files 2 and 3 are necessary to transfer. 


L13 





(zg) On NIC-80, push "STOP", depress "START" and push 
"EXECUTE" on completion of program loading. 
3. Tape Recorder 
Caer Turn ‘switch. £60 3.75 1Ps 
(b) Switch on DC bias boxes. These should be turned 
off during non-operating periods to avoid draining batteries. 
(c) Set frequency control knob on appropriate channels 
to FM and lke bandpass for signal and trigger. Voice channel 
requires 37.5 IPS and 100 Ke for operation. 
td) Turn on oscilloscope. 
4. Pulse Generator 
(a) Turn on power 
(b) Set to external trigger pulse. 
(c) Set width to 10 microseconds. On tape recorder, 


depress "DRIVE" and adjust oscilloscope for optimum display. 


B: “BXECUTLION 
1. HP-9825 
(a) Load track 0, file 8 and push "RUN" 
(b) Data entry; 
Wavelength, e.g., .6328 
OPEL eCs, e.2 7.054072 
Dwelltime, e.g., 95 
Number of scans, e.g., 125 
| Push "CONTINUE" after each entry. On last "CONTINUE" NIC-80 
will enter Display/Standby mode. HP-9825 will display 


"CONTINUE WHEN READY". 
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2. NIC-80 

(a) Verify green control light is on 

(b) Adjust CRT for optimum display of scans 

(c) Adjust DC level by knobs on bias boxes. Use low 
resolution on vertical display scale, and adjust signal to 
center around horizontal scale. 

(d) Adjust horizontal position of signal by pulse 
generator delay knob so as to center around vertical axis. 

3. To start data collection, push "CONTINUE" on HP-9825. 
When all scans are completed, NIC-80 will return to Display/ 
Standby mode and may be taken off-line. 

4. Data Analysis (HP-9825 only) 

(a) On completion of data collection, "HISTOGRAM AND 
CONTINUE" is displayed. Check that plotter and printer are 
ready before pushing ‘'CONTINUE". 

(b) Data entry; 

Date, e.g., 021577 (month, day, year) 
Time, e.g., 1423 
Cn arg, e.g., 8.07 
Power, e.g., -8 
Push ''CONTINUE' after each entry. 

(c) Unnormalized values for MTF spatial frequencies 0 
and 1 are printed for verification of absolute and relative 
values. 

(d) Data entry; "No. of class intervals" is displayed. 


Enter number of histogram bins per unit standard deviation. 
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Normally 3. Push '’CONTINUE". 
(e) Display "FREQ PLOT, YES-1". Push "1" and "CON- 
TINUE" if plot of MIF mean + 1 standard deivation for spatial 
frequency 0 - 19 is required. 
(£) Display "LN DATA? YES-1". Enter 1 and push 
"CONTINUE" if data analysis of normal distribution is completed. 
Display "'LN ANALYSIS". Push "CONTINUE" if data 
is to be converted to Ln data. 
(g) Display "DATA SUM? YES-1". Enter 1 and push 
"CONTINUE" if data summary is required. 
Display "SUMMARY, CONT.". Push "CONTINUE" to 
start execution of data summary program. 
(h) Display "SPAT FREQ". Enter spatial frequency 
number (0 - 19) for which histogram will be plotted. 
Display "REPLACE PAPER AND CONTINUE". Push 


'CONTINUE' when plotter is ready. 
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Lf X= legto “SPARD" 

maecs:par lefmt 1,£210.5:enp "wavel"”,W, “opt",0O 

enp “DWELLTIME",r0,"# OF SCANS",N 

20>+P:dim F[P],A{P],M{P],S[{P],X{P]l;wrt 2.1,W,0,r0,N;pdar 0 
for J=l to 20;rdb(2)+A;rdb(2)+B;512*B+A>+F{J] 

1f F({J)]>524287;F [J] -108576+F{J] 

SMe Xt wo EE RK L 

>: enp "CONTINUE WHEN READY",r0 

wtb 2,r0 

: for I=l to N;sfor J=1 to 20;rdb(2)+A;:rdb(2)*B;512*B+A+A[J] 
O: 1f A[J]>524287 sA[J] -10435756+A[J] 

ll: rdb(2)7*A:rdb(2)+B; 512*B8+A>r0 

P22 te bOS>524237-250- 10465754250 

13: A{J]*A[J]+r0*r0+r0;/r0+A{[J] 

14: A{J]/F[J]*A[{J] :A[J]*A[J]+S[J]+S[J] sM[J]4+A[JT]-M[J] 
Pas ext alence L,Al* | 

16: next I 

Vis ase “HISTOGRAM ?,and continue”"+sto -£lt 5 

18: eno "DATE" ,r99,"TIME",r98,"CN ARG",r94, "POWER",U 

192 “PLOE"3£Or O=1) to 20strk 1 

20: (SiOlRH4M{O]*MIO}] /N)/7/(N-=1)+S(0] -M [0] ANe4 (0) 57S(0)]=S 19) 
zi: nexe. 0 

os Ete oeort See Ome), Me M2) sit ro 7=lLegto: “pS 
23: enp “NO CLASS INTERV", T;beep 

242 START" -enp "FREO PLOT, YES=L" ,r9 2712 r92=1l-9to ~FREPLO” 
o>: Ir E97 —-i gto fps" 

26; “EV tenp “EN DATA, YESa=s, co7eut r97=lrleZegto Fin" 
2/3 "DS s0+kreno “DATA SUM? YES-l"oAs16 KA=letrk O¢ldrt 9 
28: "“PLOTO":senp “SPAT FREQ",A;A+1l*-OQ;trk 1 

29: asp "REPLACE PLOT PAPER AND CONTINUE";been;stp 

S02 (Ol s3"55 10)-6r (Ol +375 -e) >D 

31: S[(Q] /T*R:C-2,2*Rer9lscsiz 1.6,1,5,l 

B2: SCIVCG 2787), “2271.70 t DyOpie etc. 10,0527 <0 aL 
B32 Plt. rode 2 0s le pt Teor 257 pe bor 1=0 LO a7 OY ai 
$4: teenie, 5, '- lol Grspee r9lFR/2=8/10;1,1 

$9: DeGaeo ek 727 10;,1,2 next 1 

363: —-3*T+r88;£xd 0 

os LOrel-Gvcosbeovy Riolt 17,01, 701. 1,-.01, 2 

38: ltr I-R/3,-.08:lbl r388 

B95 EOOt roses mexe Ioan esiz 1; 7S.) 

AQ; Ltr -Ce3*k, te6- lol "RELATIVE DISTRIBUTION OF <“2r™ 


PO mnAID UO & Wh fF © 


Figure 8A. Data Analysis Program Listing 
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pen; ple M(O],Q9;plt M[Q],.7,2 

ltr CHR, = sls lbl “UNETS = .1/",T,™ STANDARD DEVIATIONS" 
"GURVE sfeor 1=C to D by “(b-€) 7/1003 (I-Mio)j7slolerl 
rl*rl/2+rl;exp(-rl)*R/(S[Q]*y¥ (2*7) )+r2 


© plt fpeZ2,-2emnext Peeer 2; rZ—-Rer3 
¢ Cer2:r2+Rer3;f£lt 5S;AeA 


for K=ll to 114+T*6:07*rK:next K 

"HIST" sfor I=l to N;Cer2;r2+R-r3 

ld£ 1,A[{*];for K=ll to 11+T*6 

fit S:1f A[Q])>r2:1if AlO]<r3;rkKt+l>rkK 
r3er2;r3+Rer3;next K;next I 

CerZe-r2tkee3:sfor K=ll to LItT*6:rkKer1 

DuyNee pens Plt £Z,0,—-2-Olt £2,r1,2,0.t £3,r-1,2 

DIG e oa0, 2 
(azaMUOW)7S(Olerisr7*ri/2>r7- exp (=27)*R7 (S10) *yalem) jeer 
(ReMi) /S Ole ier tei 727 exp (—r7) *R/(Si0)]* 4 (2*0)) 465 
R*r4t+abs (r5—-r4)/2*Rer3;:r8/Rers 

(N¥*r8—-r1l*N)7 2/(N*r8)-r93r94+r1lU+rld 

r3er2;r3+R*r3;next K 

“LABEL” sC+R/ 3+rl-£xd O:ltr rl,1.45 

lbl "SPATIAL FREQUENCY",Q-1 

if r97=l;ltr rl,1.25;lbl "TEST DIST LN NORMAL";gto +2 
Ltr rl,l,Z5) fol “TES! ois) VoOnwAL,” 

Lepr ee ler. 35. tol “NO OF SGANs "nN 

Exd 4-1tr el 1.15? 1b] “WAVELENGTH”, a: ltr rl1,1,05 

loi OR TTeS” 0 

EXGeO tere tl, «9 lb) “BATE "139, * TEMG .~r9serilec 5 


: MOIR] S=eleGr rifle sos tek "St, dev",;S10] 


Ltr rl, 1.453151 “MEAN =o] 

Ler iit Zoo). “CHE SO", 7 lO; ltr ri, i.sl>->ixd 3 
lbl 7G". r94:txd 0 

Lo Ve au 

for I=l to 90;0+rI;next I 

LE L9/=legto “DS” 

gto ss) ins 

"IN ":dsp "LN ANALYSIS";stp ;trk 1 

ror J=l to 20;0+M[J]+S[J] ;snext J 

for I=l to N 

IGE T,A1*);for J=l to 20e1f A{J]=071£F J=20;gto +5 
1f A{[J]=O:next J 


Figure 8B. Data Analysis Program Listing 
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SSS 


In(A[J])+AT[J] 


>: ALJ]*A[J]+S[J]+S[J] sM[J]+A[J] +M[J] 


next J:rcf I1,A[*] 

next beagto = PLOT" 

oR ePEOewce | —3420.,— 12 ome. oseeed. 0 -eerz 1615, 1 

DEC 20,0 OL Upc oo, oie Sho oe ks. 2 

fOr t= tO okt l-  wUlehe ort I—-l,—. 0c nexcte 1 

LOL t= tCOsZemoy S5Sslsleros+ ltr 195-13, -s059+tbpl 695 
next, f-ixda 1 

EOnmel—Umee 1 DY. ber —]-65, 121 bl te ple--). 9,172 

EOr 1=0 to Jeaby .heetr —l3>, 1) dl Te olt —-1..9,1, 1 

Dit. —1.,6).b9¢2- meee ot 

Ltr 2,-~-b9- 1b “SPATIAL FREQUENCY IN L/ MILEIRADIANS" 
ltr 2,1,.45;lbl "MTF MEAN AND + —- 1 STANDARD DEVIATION" 
fxd ouster 24....35 

ble "“DATE",r99," TIME"..r98," WO OF SCANS”™..4 

EExd 427 .tr 5,).25-¢1D “WAVELENGTH ',N,. “OPTIEGS”,.0 
ERGCwoe Ltr 5, .sL>efed Sylolt “Gy, roger sd OFfbl "EY, 


: M[1l]+r93 


Corey, IU. Ler S21 peo Ol UME 

for Q=1 to 20;50+90er43 570+0>r49 sM[0] fr 93-743 
5[Q] /r93+r49 

pit OS) hpas- E49, lsolt Cal peastr4?, 2 


s OLt. O-14057646,. Olt 0O-.75,045,2 Den next O-gto. "ft ” 


Figure 8C. Data Analysis Program Listing 
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Mm OONTAIU BW MNF © 
ee ee ee 


© ae ee ee 
ee 


tis 


Lo 
a) 
ee 


dsp “SUMMARY ,CONT";beep;stp 

"X"sftor Q=l to 20 
S10) 7828 
M@AwSs co. eCs MIO) +3 *5(0\+D 
for K=ll to 11+T*6;0-rK:next K 

"HIST" :sfor I=l to N:Cer2;:r2+Rer3 
=) eee 
fats1,A|*\|-tor K=l!1 to Lit+tT*s 
Eiteo-e dt AlOl> E21 Ora rktlerk 
r3er2;r3+R*r3;next K;next I 
Cer2:r2+Rer3;for K=ll to L1lt+T*o;rkerl 
rl/Nerl 
(r2—M[0)])/S(Qler7sr7*r7/2-r7;exn(-r7)*R/(S 

(esa Ole le titans oy ec) gexe(-r7)*R/(S(O]*¥ (2"0))*r> 
R*r4+abds(r5—-r1r4)/2*Rerd;3 r8/Rer8 

(N*r8—rlL¥*¥N)° 2/(N*r8)%r93r94+rl0-rld 

cr3er2:r3+Rer3snext K;rlO+A [QO] :0erlO;next OQ;gto "PRINT" 
VU Em ese ok, ae eX FOO, ox ys oh ie eee oe OC 

Pit Soke eM ee xt mee seek, Ft Oo eek, 6 CLC S [2k 78O,4 72s 
Emme 3,455, "DATA SUMMARY" 

tise taaeO "L1" 
WE, O 65,5 

web 5,170,190 

wet O,.2,099,r98 
Web yo. U0l 

WE tsO 4b pr 9d, , O 
Web 6,20, 10 

Piet OX, TES. DIST .NORIAL” ,25x. "TEST Dis? LN NORMAL” 

wrt 0,4 

Wie Or kO 10 

fic Seeeoe ee oe AL ER. ox p “MEANY Aix "SYD EVY 7 20x, GCHiL oO 710s, “Gal 
Wet 6 to web, LU, LO 

iG) O 7, ox eo, Ose lL. 55x el... 5, ox VELL. 

tor I=) tO .20-wre 6 76,..—17411),S(1],4)1lenext I 

FOr [=] tO 20-Wtb o,27,10-next [ 

Peel go 2xze ll ss 

if Zell gto "Ll" 

tre Oar 3 

Molermte fyi 2x ,el 5 

EOGml=teto 2ZO* wre 6.7 -4 Lt) next 1 

Erie o lar € 


Figure 9. Data Summary Program Listing 


120 





110600 
4304 
44034 
10377 
492340 
0 
2175677 
4304 
3001737 
3001745 
1916 
3001737 
3001745 
1614 
Z00L737 
3001745 
1600 
3001737 
3001745 
1620 
110020 
2405677 
2001500 
STENS 76 
10777 
6474 

et 

4473 

31 ELS 6 
403031 
6474 

a6 

4473 
2125576 
2703677 
27 
3001742 
1600 
3001715 


Figure lOA. 


MEMA(600 /READOUT LIGHTS 
LCWORD /CONTROL WORD 

STU /READ STATUS 

ANDA(377 /COMPARE SETTINGS 
A-MZ(340 /CHECK HARD SETTING 
JMP Q /RETURN CHECK LOOP 
ZERMA 1677 =/0 TO ACC & 1677 
LCW(O /LOAD CONTROL WORD 
JMS@INP /INPUT 

JMS@1745 /PUTAC 

1616 /ADR OF WAVELENGTH 
JMS@INP /INPUT 

JMS@1745 /PUTAC 

1614 /ADR OF OPTICS 

JMS@INP /INPUT 

JMS@1745 /PUTAC 

1600 /ADR OF DWELLTIME 
JMS@INP /INPUT 

JMS@1745 /PUTAC 

1620 /ADR NO OCF SCANS 
MEMA(20 /COUNT VALUE 

ACCM 1677 =/COUNT INDEX 

JMS 1500 /CAL ADRESS 
MEMA@1576 JCA VALUE 
ANDA(777 /9 BITS 

HSP /IS HP 9825 READY 

IMP 3L /WAIT LOOP 

PHS /WRITE TO HP 9825 
MEMA@1576 /CAL VAL RECALLED 
RIS lL /SRIGHT SHIFT 

HSP /IS HP 9825 READY 

JMP 36 /WAIT LOOP 

PHS /WRITE REMAINING DATA 
MPOM 1576 /INCR, CAL ADRESS 
MMOZ 1677 =/DECR.&TEST 

JMP 27 /RETURN LOOP 


JMS@1742 /GETAC 
1600 /ADR,. OF DWELLTIME 
JMS@1715 /FIX 


Nic 80 Program Listing 
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3.1422 MEMA@1722 /FACM 
4311 LOWELL /LOAD DWELL TIME 
110000 MEMA(O /0 TO ACC 

4312 LDELAY /LOAD DELAY TIME 
2A ie MEMA 1711 /STAND BY LOOP 

4304 LCWORD /LOAD CONTROL WORD 

4302 SET /SET MEASURE FF 

S220) STOP /READ 1 SCAN 

6464 HSRF /SKIP ON HP REQUEST 

53 JMP 53 /RETURN STAND BY LOOP 
44463 RHSR /READ HP 

462314 A-MZ (314 /L FOR FUTURE BRANCH 
64 JAP 64 /FUTURE BR.ADR, 

3001723 JMS@1723 /SHORT TERM RESET 

2111704 MEMA 1704 =/B4(0) 

2405677 ACCM 1677 / 

2111626 MEMA 1626 /# OF SCANS 

2405672 ACCM 1672 / 

PAE a MEMA 1712 /CONTROL WORD #2 

4304 LCWORD /LOAD CONT.WORD 

4302 SETM /SET MEASURE FF 

52.20 STOP /READ 1 SCAN 

100 JMP 100 / 

0 Q /NOT USED 

Q Q /NOT USED 
110200 AEMA( 200 / 

4304 LCWORD /DISPLAY CONTROL 
2111700 MEMAL L700 /S80%0), BEET £ND 
2405677 ACCM 1677 / 
27h 2701 MMOA 1701 /BO(1777),RIGHT END 
2405076 AGEM 1676 7 

110040 MEMA(40 /32 POINTS 

24056075 ACGY 16755 / 

25674 ZERM 1674 / 

Seo 7 MEMA@1lO77 =/INTEGRAL BASE LINE 
2505674 A+MM 1674 / 

3L1L1676 MEMA@1676 / 

2505674 A+MM 1074 / 

Figure 10B. Nic 80 Program Listing 
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Zio 
2705676 
2707675 
ilk 
2111674 
5026 
2405674 
ZEEE O9 
2405677 
3 G7 7 
2471674 
3405677 
ZEST 7 
2463701 
25 
PLO Z 
S00R722 
2111703 
2405677 
2111700 
2405676 
3165676 
3165677 
2125677 
2135676 
2463701 
142 
3001776 
3001622 
4032 
2703626 
71 

23 


Figure lOC. 


MPOM 1677 / 
MMOM 1676 / 
MMOMZ 1675 / 
JMe i) 7 

MEMA 1674 / 


RASH6 /DIVIDE BY 6 

ACCM 1674 /AVERAGE DC COMPONENT 
MEMA 1700 /BO(Q) 

BCCM lows 7 

MEMA@1677 / 

A-MA 1674 /SUBTRACTING 


ACCM@1677 / 

MPOMA 1677 / 

A-MZ 1701 / 

IMP 126 /RETURN ADR,LOOP 
MEMA (2 
JMS@1722 
MEMA 1703 
ACCM 1677 
MEMA 1700 /BO 
ACCM 1676 
ZERM@1676 
ZERM@1677 
MPOM 1677 
MPOMA 1676 / 
A-MZ 1701 / 


/MOVE BO+B2 
/B3 ZERO LOOP 


~ 


eS 


JMP 142 /RETURN ADR, LOOP 
JMS@1776 / 
JMS@1622 /DFT DATA TRANSFER 


ASCNTR /ADVANCE SWEEP COUNTER 
MMOZ 1626 /DECREASE SCAN CCUNT 
JHP 71 /RETURN ADR,NEXT SCAN 

JMP 53 /RETURN ADR.STAND BY L, 


Nic 80 Program Listing 
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E52 
110024 
2405675 
211702 
2405677 
Ze 3 
2405676 
STErOd 7 
10777 
0474 
1456 
4473 
SLEE677 
405031 
6474 
1463 
4473 

3 Eos 
LOdd 7 
6474 
1470 
4473 

3 Ae 6 
405031 
6474 
1475 
4473 

Zi 2507 1 
2125676 


1001445 


Figure LOD. 


152 /RETURN ADRESS 
MEMA(24 /# OF SPAT.FR,. 
ACCM 1675 /CONTING INDEX 
MEMA 1702 /B2,REALS 

ACCM 1677 /CURRENT B2 
MEMA 1703 /B3,IMAGINARIES 
ACCM 1676 /CURRENT B3 
MEMA@1677 /VALUE B2 
ANDA(777. /9 BITS 

HSP /HP READY 

JAP 1456 /WAIT LOOP 

PHS /WRITE TO HP 
MEMA@1677 /RECALL B2 

RIS ll /RIGHT SHIFT 

HSP / 

IMP 1463 / 

PHS / 

MEMA@1676 / 

ANDACI77  / 

HSP / 

IMP 1470 / 

PHS / 

MEMA@1676 / 

RIS ll / 

HSP / 

JMP 1475 / 

PHS / 

MPOM 1677 /INCREASE 32 
MPOM 1676 /INCREASE B33 
MMOMZ 1675  /DECREASE & TEST 
1454 /RETURN ADR. 
JMP91445 /RETURN TO MAIN PR, 


Nic 80 Program Listing 
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